Работа посвящена разработке, реализации и исследованию алгоритма расчета частотно-зависимых времен пробега сейсмических волн. Алгоритм состоит из двух модулей - сглаживания и трассирования лучей. Сглаживание производится по зоне Френеля и зависит от положения источника, частоты и поля времен. При трассировании лучей используются производные времен пробега, полученные методом быстрой маршировки для решения уравнения эйконала. Разработанный алгоритм исследуется на реалистичной синтетической модели, построенной по реальным данным из Юрубчено-Тохомского месторождения. Для этой модели проводится сопоставление результатов работы разработанного алгоритма, лучевого метода и волновых полей, полученных конечно-разностным моделированием.